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We review the out-of-equilibrium properties of a self-gravitating gas of particles in the 
r~| ' presence of a strong friction and a random force (canonical gas). We assume a bare 

^ ' diffusion coefficient of the form D{p) = Tp^/" , where p is the local particle density, so 

1^ ' that the equation of state is P{p) = D(p)p. Depending on the spatial dimension d, the 

G , index n, the temperature T, and whether the system is confined to a finite box or not, 

" ■ the system can reach an equilibrium state, collapse or evaporate. This article focuses on 

the latter cases, presenting a complete dynamical phase diagram of the system. 

ly^ Keywords: Gravitational collapse; stochastic processes; dynamical phase transition. 

1. Introduction and background 

^ I There is currently a renewed interest in the study of long-range interacting systems, 

^ . particularly outside the realm of astrophysics. These systems may exhibit inequiva- 

Q I lence of statistical ensembles {e.g. canonical vs microcanonical) which affects their 

equilibrium properties (possibility of negative specific heat in the microcanonical 
ensemble) and even more dramatically, their dynamical properties. In the context 
of self-gravitating systems, the Newtonian dynamics is too complicated to permit an 
exhaustive analytical treatment. Hence, in the past few years, we have developed a 
model of self-gravitating particles in the presence of a strong friction and a random 
force, ^ for which inertial effects are negligible. The main interest of this model is 
to be analytically tractable in many situations, while presenting dynamical phases 
I reminiscent of their Newtonian counterparts. In addition, this model is intimately 

related to the Keller-Segel model of bacterial chemotaxis.^ 

We thus consider particles obeying the equations of motion ^ = —\7^+\^2Dr]i, 
where iji is a delta correlated random Gaussian force, and $ is the gravitational 
potential. In a proper mean-field limit, which becomes exact for an infinite number 
of particles, the density obeys a Fokker-Planck (or Smoluchowski) equation coupled 
' to the Poisson equation (the gravitational constant is set equal to G = 1): 

^ =V-(VP + pV$), A$-p. (1) 
ot 

The pressure is related to the diffusion coefficient by P{p) — D{p)p, and the isother- 
mal case corresponds to D = T, while we will consider here the more general poly- 
tropic case D = Tp^^"^.^ In addition, the system of total mass M can be placed in 
a spherical bounded domain of radius R or in an unbounded space. 

In the context of chemotaxis,^ p is the density of bacteria, and c = — $ is the 
density of a chemical that they secrete. The bacteria are attracted by the regions 
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of high density of this chemical which generates an effective long-range interaction 
between them. This interaction exactly takes the form of gravity when neglecting 
the diffusion and the degradation of the chemical. 

Although wc shall sec that the actual phase diagram depends crucially on the 
value of the index n, the general physics is controlled by the parameter T which we 
assimilate to the temperature. When T is small, the kinetic pressure is not strong 
enough to compensate gravity and the system may collapse. For large T, the system 
is at equilibrium in a bounded domain or evaporates in an unbounded domain. In 
the next section, we present the main results that we have obtained in^"* concerning 
these different dynamical phases. 

2. Collapse and post-collapse dynamics and the critical index n* 

Defining the critical index n* = -^^12' system can be shown to always have a 
polytropic equilibrium state for < n < n*, whether the system is confined or not. 
Hence we shall concentrate on dynamical properties of the case n > n* and n < 
(including the isothermal case n = oo (d > 2)), and on the interesting particular 
case n = n* (including the isothermal case n = oo (d = 2)). 

2.1. n > n* and the isothermal case n = oo (d > 2) 

We have shown-'^'"^"^ that in a bounded domain, and for T < T^, the system displays 
a finite time singularity at which the central density diverges as po{t) ~ (icoii — t)~^. 
The density profile takes a scale-invariant form p{r,t) = po{t)f[r/rQ{t)\, where the 
core radius r^(t) is related to the central density by ro(t) = T / p\ ^^"(t). The density 
scaling function decays as f{x) ^ x^", with a = '^'^i s-nd can be analytically 
obtained in the case n = oo (d > 2).^ In this important case, <con was also computed 
exactly close to Tc.^ 

For t > tcoU, during what wc dubbed the post-collapse regime!,* a delta peak 
condensate grows at the center as Mq ~ (t — tcoii)'^^""^^^'-^"'"^, and ultimately 
saturates exponentially to the total mass M, with a rate which was computed 
analytically for low T, by using semi-classical quantum methods.'' Meanwhile, the 
residual density obeys a reverse scaling of the form Prcs.(^iO = Pr{t)fr[r/rr{t)], 
where Pr{t) ~ (i — tcoii)~^i and fr decays with the same exponent a = as /. 

In an unbounded domain,^ the gas collapses below a non universal critical tem- 
perature Tc depending on the initial state, or evaporates for T > Tc. In this evapora- 
tion regime, gravity becomes asymptotically negligible and free diffusion is observed, 
with sub-corrections due to gravity which have been obtained analytically.^'^ 

2.2. n < 

The case of negative n.^ including the logotropic case n = — 1, is similar to the case 
treated above for n < —d/2, with a collapse scaling exponent a = However, 
for —d/2 < n < 0, the collapse is controlled by the T = fixed point^'^'^ with 

" d+2- 
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2.3. n = and the isothermal case n = oo (d — 2) 

In this critical case,^'^'^ Tc is independent of the radius R of the confining box 
(Tc — M/4, in d = 2). For T < Tc, there is still a finite time singularity but a 
delta peak condensate of mass Mc — {T/Tc^/^M forms at the center, at t = icon- 
Moreover, the residual density takes a pseudo scale-invariant form with an effective 
scaling exponent a{t) very slowly saturating to a = d, a.s t goes to tcoii- In an 
unbounded domain, the gas collapses below the universal critical temperature Tc, 
or evaporates for T > Tc- In this evaporation regime, the scaling density profile is 
strongly affected by gravity and the gas expansion is such that {r'^{t)) ^ {T—Tc)t.^''^ 

3. Summary 

In TablelU we summarize the static and dynamic phase diagram of a self-gravitating 
gas of particles with a bare diffusion coefficient D = Tp^/^^, where T is the temper- 
ature and d the spatial dimension. This table illustrates the crucial role played by 
the critical index n* = 



Table 1. Static and dynamic phase diagram of a self-gravitating gas. 



Index 7\ 


Temperature 


Bounded domain 


Unbounded domain 


n = oo (d > 2) 


T > Tc 


Mctastablc equilibrium state 
(local minimum of free energy): 
box-confined isothermal sphere 


• Evaporation : 
asymptotically free normal 
diffusion (gravity negligible) 
• Collapse: 
pre-collapse and post-collapse 
as in a bounded domain 


T <Tc 


Self-similar collapse with a = 2 
and self-similar post-collapse leading 
to a Dirac peak of mass M 


< n < Tl* 


T > Tc 


Equilibrium state: 
box-confined (incomplete) polytrope 


Equilibrium state: 
complete polytrope 
(compact support ) 


T < Tc 


Equilibrium state: 
complete polytrope (compact support) 


■n^ < n < <x, 
11 < 


T > Tc 


Metastable equilibrium state 
(local minimum of troo energy); 
box-confined polytropic sphere 


• Evaporation: 
asymptotically free anomalous 
diffusion (gravity negligible) 
• Collapse: 
pre-collapse and post-collapse 
as in a bounded domain 


T < Tc 


b!clf-slm!lar collapse with a = -^n/in ~ 1) 
or Q = 2d/(d -1- 2) (-d/2 < n. < 0) and 
post-collapse leading to a Dirac peak 




T > Tc 


Equilibrium state: 
box-confined (incomplete) polytrope 


Self-similar evaporation 
modified by self-gravity 


T < Tc 


Pseudo self-similar collapse 
loading to a Dirac peak of 
mass (T/Tc)^^^ M + halo. 
This is followed by a post-collapse 
leading to a Dirac peak of mass M 


Collapse 
as in a bounded domain 


T = Tc 


Infinite family of steady states 


As in a bounded domain 
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